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trans-GZK Cosmic Rays: Strings, Black Holes, Neutrinos, or all Three? 
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We review the scenario in which "strongly interacting neutrinos" are responsible for inducing airshowers with 
inferred energies E > 8 x 10^® eV. This possibility arises naturally in string excitation models having a unifi- 
cation scale effectively decoupled from the Planck scale. We then show that phenomenological quantum gravity 
considerations reveal an equivalency of "mini-black hole" and strongly interacting neutrino pictures for explain- 
ing trans-GZK events. This equivalence can be exploited to predict single particle inclusive distributions. The 
resulting observable consequences in airshower development are studied using the Adaptive Longitudinal Profile 
Shower (ALPS) simulation. 
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1. INTRODUCTION 

The publication of the seminal papers by 
Greisen and Zatsepin & Kuzmin P predicting 
the energy attenuation of protons overs tens of 
megaparsecs due to interactions with CMB pho- 
tons (GZK effect), together with the observation 
of cosmic ray airshowers with inferred energies 
of > 10^° eV, has created a vigorously pursued 
area of particle astrophysics. In fact, as of 2004 
there are many more papers proposing explana- 
tions for the existence of these trans-GZK air- 
showers than there are recorded events! Due 
to this relatively sparse sample of events, open 
questions include the existence of statistically sig- 
nificant clustering/anisotropy, correlations with 
known astrophysical source distributions (e.g., 
QSOs), the composition/charge of the primaries, 
and whether the production mechanism is top- 
down or bottom-up. The present paper, delivered 
at the 2004 Cosmic Ray International Seminar 
(CRIS 2004), summarizes a bottom- up scenario 
whereby the primary incident on the Earth's 
atmosphere is a neutral, non-hadronic particle 
such as a neutrino that generates a hadron-like 
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airshower due to the primary interaction being 
above a low-scale unification threshold. Further 
details concerning the theoretical considerations 
sketched here, additional figures and tables, and 
more complete bibliographies can be found in the 
referenced papers. 

2. THEORETICAL MOTIVATIONS 

If trans-GZK events originate from sources far- 
ther away than 50 — 100 Mpc, the most likely 
candidate for the primary is some type of neutral 
particle. A Standard Model neutrino can propa- 
gate over cosmological distances with little energy 
loss, but, of course, the interaction cross-section 
<JuN~>tx is several orders of magnitude too low to 
generate the observed showers. However, it has 
been conjectured for some time that perhaps a 
"strongly interacting neutrino" could exist due to 
some as yet undiscovered "new physics" , and thus 
possess the required interaction strength at the 
appropriate CM energies With the realization 
that higher-dimensional string theories may allow 
interaction unification at energies many orders 
below the standard 4-dimensional Planck scale of 
~ lO-'^^GeV P], interest in the strongly interact- 
ing neutrino picture has been revived, and has re- 
sulted in new work involving rigorous theoretical 
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considerations within the framework of a specific 
model 01 . Additional features relative to cosmic 
ray airshowers have been extracted, although on- 
going work continues to refine the predictions rel- 
ative to conventional airshower observables 5 . 

2.1. Revievif of the String-inspired Scenario 

More specifically, if the CM energy of the 
neutrino-nucleon is above a "low energy" uni- 
fication scale, the neutrino-quark interaction is 
"strong" , and a leptoquark resonant state can 
form, a,yN^LQ-[os^ix- In establishing model 
properties using the phenomenology of higher- 
dimensional string theories, it is worth noting 
that "typical" tree level amplitudes for strings 
yield a cross section that is too small to ac- 
count for the trans-GZK events 0. This does 
not rule out the present model, but does empha- 
size that weakly-coupled string theories are in- 
adequate here, and one must calculate with the 
strongly-coupled theory (and non-perturbatively 
when possible). This is analogous to attempting 
to calculate quark interactions using QCD. 

The basic building blocks of the model are 
unitarity of the S-matrix, a rapidly (exponen- 
tially) rising level density of resonances, unifi- 
cation of interactions with strength ~ Standard 
Model strong interaction at string scale M*, and 
duality of resonances in a given channel with 
Regge exchanges in crossed channels. Express the 
total cross section in terms of partial waves and 
absorption coefficients using the optical theorem, 

Ovq->LQrcs{s) = ^ (2j + 1)(1 - 'qj C0s(2(5j)). 

3 

Note that as s increases, the resonances are no 
longer purely elastic with absorption slowing the 
growth of the total cross section, and where 
N{s) is the resonance level and equals the max- 
imum angular momentum; ignoring corrections 
0{Mz /M^,), avq-^LQrcs tends to a constant value. 
The level density of resonances is well represented 
for the first few resonances by 

dNo oc l.2mQ « exp l.2A{so/Ml) . 

2.2. Strings and Quantum Black Holes 

A connection between string theories and quan- 
tum black holes is not a new idea, and there is no- 



table recent work in this area T . Extending this 
to the case of interest here, Domokos & Kovesi- 
Domokos have demonstrated a non-perturbative 
equivalence of the string excitation model with 
that of quantum "mini black holes" 

The treatment in reference ;8 is based on a 
statistical mechanical analysis of the previously 
introduced string model beginning with consider- 
ation of the microcanonical density matrix of the 
final state, 

p = J2\N,a}{N,a\S{e~ NM,) . 

a 

A straightforward manipulation yields the stan- 
dard passge to the canonical ensemble with well- 
defined entropy 5' and a temperature that asymp- 
totically approaches the Hagedorn temperature, 
Th — M*/37r. Recall that in modern QCD, the 
Hagedorn temperature Th is interpreted as mark- 
ing the deconfining phase transition from the low 
temperature hadron phase (quark confinement) 
to the quark-gluon plasma. In terms of modern 
quantum gravity and string theory, the entropy 
of a d-dimensional quantum black hole (QBH) 
equals the string entropy S at Tqbh ~ ~ Th- 
It is straightforward to derive additional relations 
between 5, M*, Th, and Rschwarz- 

The asymptotic estimate then has important 
consequences in allowing a non-perturbative cal- 
culation of the single particle inclusive distribu- 
tion similar to the statistical mechanical analy- 
sis carried out at the tree level for specific string 
models by Amati & Russo (1999). Naturally, cal- 
culated quantities are a function of the number of 
"extra" dimensions in the model (i.e., the num- 
ber of dimensions in addition to the conventional 
4-dim space-time). Depending on model choices 
for characteristic scales, the resulting multiplici- 
ties can be similar to those expected in reactions 
initiated by a heavy nucleus 0. 

2.3. The Transition from SM to Unified 
Regimes 

Because the string-QBH model must contain 
resonances with the number of states an expo- 
nentially growing function of the resonance mass, 
there is a (nearly) 6'-fcn transition from the SM 
to the unified regime. However, there is obviously 
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Figure 1. Smoothing function used for transition 
to unified interactions regime showing the nearly 
step-function behavior due to the exponentially- 
rising level density. 



Figure 2. Integration of the approximately step- 
function neutrino-quark cross section over the ap- 
propriate PDFs; the dashed lines indicate the ap- 
proximate region that would generate observable 
deep showers. 



no exact 0-fcn transition in Nature, thus imply- 
ing that there must exist a finite width transition 
region. With no analytic solution describing this 
transition available from the physics, we choose 
to treat it empirically by adopting a convenient 
mathematical smoothing function to represent 
the neutrino quark cross section. Nevertheless, 
the chosen form exhibits the desirable physical 
characteristics of nearly step-function behavior 
saturating at a TBD strong interaction strength 
while remaining continuous and without violating 
unitarity, 
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In the above, the conversion from the Lab frame 
to the CM frame is via s = 2mNucE, and 
A^o = 2mNucEo/M^ such that at E — Eq we 
have a^q^LQRos{E = £'o) = \a^q{saturated). 
The choices of C and Nq determine the value of 
duq{saturated)^ taken here to be 50 — 100 mb. 
A typical example is shown in Figurenfor a string 
scale 70 TeV. Of course, the neutrino-quark 

cross section must be integrated over the parton 
distribution functions (PDFs) of the participat- 



ing partons in the nucleon to arrive at the final 
(observed) quark-nucleon cross section, 



(^vq^LQT:ca^ix{s, M^,Nq) 



dxf{x)a{s), 



for s = xs and momentum fraction x. The in- 
tegration is over valence quarks only since there 
are no leptogluons in current string models, and 
the sea quarks contribute primarily around x = 
0. Figure |21 shows the result of integrating the 
a^q^LQ.ea of Figure n using the CTEQ6 PDFs. 
There is a region where the cross section would 
give rise to deep showers, and this is indicated 
(approximately) by the dashed lines in the figure 
(this is further discussed below). 

3. OBSERVABLE CONSEQUENCES 

Thus, as developed in the previous Section, the 
strength and properties of the unified interaction 
can lead to a hadron-like airshower initiated by 
a specific neutral string/QBH state (here taken 
to be a neutrino), but with some potentially ob- 
servable differences compared to those generated 
by hadrons or nuclei. We study the generated 
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Figure 3. Comparison of X^ax values generated 
by ALPS, CORSIKA-QGSjet, and the Bartol- 
QGSjet simulations for proton primaries. The 
error bars shown are from ALPS. 
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Figure 4. ALPS simulations of 1000 airshowers 
generated by 3 x 10^° eV protons and (strongly 
interacting) neutrinos. 



airshower properties using the Adaptive Longi- 
tudinal Profile Shower (ALPS) simulation origi- 
nally created by P. Mikulski ALPS is similar 
to other hybrid simulations that use subshower 
parameterization instead of tracking every indi- 
vidual particle. Its performance is comparable 
to CORSIKA-QGSjet and the Bartol-QGSjet Hy- 
brid Simulation, as shown in Figure 01 for average 
Xmax values generated by proton primaries as a 
function of primary energy. Similar results hold 
for Nmax values. 

Using the same string model parameters lead- 
ing to Figures ^ and |21 Figure 01 presents ALPS 
results showing that, in general, a proton will pro- 
duce more electrons in a given airshower than a 
strongly interacting neutrino of the same energy 
although the depth at maximum electron produc- 
tion is quite similar. The reason for the Nmax 
difference is that a larger number of prompt lep- 
tons (mostly muons) are produced by the neu- 
trino. However, from Figure [51 it is also seen 
that a 4.5 X 10^° eV neutrino has both the ap- 
proximately same Xmax and N^ax distributions 
as the 3 x 10^° eV proton, so that it would be vir- 
tually impossible to distinguish airshowers gener- 
ated by one or the other. In these cases, an addi- 
tional airshower observable such as muon number 




Figure 5. ALPS simulations of 1000 airshowers 
generated by 3 x 10^° eV protons and (strongly 
interacting) neutrinos and 4.5 x 10^° eV neutri- 
nos. The higher energy neutrino events ('-I-' sym- 
bols) populate essentially the same region as the 
proton events ('o' symbols), and are not easily 
distinguished in this plot. 
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Figure 6. Preliminary extrapolations from ALPS 
simulations for the variation of X^ax and Xq in 
the SM unified transition region for a 70 TeV 
string scale. 



4. DISCUSSION AND SUMMARY 

It should be mentioned that the idea that 
mini black holes might be produced in extra- 
dimensional theories with energy scales orders of 
magnitude lower than the Planck scale has been 
explored previously by several authors (lOj . How- 
ever, the model presented here differs in several 
important aspects: 

1. "TeV-Scale Gravity is set by the SM elec- 
troweak scale (EW), 

2. Interactions are not unified in such models 
(not a GUT), so there are no strong scale inter- 
actions, 

3. The cross sections for the black hole in- 
teractions are derived using a semi-classical ap- 
proach utilizing geometric total cross sections, 
aBH ^ ''^^'schwarz' ^^'^ resouauccs are not in- 
cluded; while these cross sections are considerably 
enhanced compared to SM EW values, they are 
still significantly smaller compared to SM strong 
cross sections. 



is required. Still, since it is expected that sources 
possess a non-monoenergetic injection spectrum, 
over a large number of observed events with dif- 
ferent primary energies, the best discriminator 
may be identification of a specific source candi- 
date since the source distance can eliminate pro- 
tons as candidate primaries. 

Another possible discrimination technique is 
the observation of highly inclined or nearly hori- 
zontal airshowers since the expected width of the 
SM — > unified transition region should produce 
an excess of observed showers relative to that ex- 
pected from protons alone. Predictions (extrap- 
olations) for how Xmax and Xq scale with the 
changing cross section in the transition region are 
shown in Figure Also note that the fluctua- 
tions around the mean values shown in Figure El 
also increase very rapidly with decreasing cross 
section. The interesting point here is that an ob- 
served flux of inclined showers can be used to di- 
rectly constrain the string scale, and, depending 
on the model favored by Nature, this may be the 
only way to do this in the next decade (or be- 
yond). 



Gurrently, experimental results do not rule out 
the model discussed here, TeV-scale gravity mod- 
els, or many other competitors. However, the 
next generation detectors, especially Auger and 
ASHRA, that are nearing completion are quite 
capable of discrminating among the possibilities. 

It is also interesting to note that a detailed 
analysis of the AGASA data appears to indicate 
that the highest energy cosmic rays with E > 
8 X 10^^ eV may be distributed on the sky differ- 
ently than those having energies 4 < < 8 x 10^^ 
eV, and this may support the idea that the high- 
est energy airshower events are generated by neu- 
tral particles that have propagated distances in 
excess of 100 Mpc 11 1. The main difhculty with 
such a scenario is that, at this time, it is difficult 
to understand how a > lO'^'^ eV neutrino is pro- 
duced without first generating a ^ 10^^ — 10^^ eV 
hadron, and it is simply not known if such an as- 
trophysical engine exists. However, it is probably 
too soon to say with absolute certainty that such 
sources cannot exist. 

Finally, we summarize the main points pre- 
sented here. Our work is continuing in both the 
theoretical and the simulation areas, and will be 
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reported in future papers. The robust features of 
the string/QBH phenomenology developed up to 
now include the following: 

1. For given values of C and A^o, the width of 
the SM — > unified transition region broadens as 
the characteristic string scale increases, in- 
dependent of the choice of smoothing function; 
a lower limit on the flux of deep showers corre- 
sponds to an upper limit on A/*, 

2. The energy at which the rising cross sec- 
tion reaches 1/2 the saturated unified value de- 
termines the presence (or lack) of a dip at high 
energies in the CR spectrum, 

3. Exploiting the String-QBH equivalence pro- 
vides a non-perturbative method for calculating 
single particle inclusive distributions, 

4. The presence of extra dimensions affects 
the single particle inclusive distributions in a way 
that is potentially observable in cosmic ray air- 
showers, and may provide a mechanism for deter- 
mining the number of extra dimensions preferred 
by Nature. 

The model also yields a natural way to include 
an additional component in the EHECR distribu- 
tion if protons, pions, and neutrinos are produced 
in astrophysical engines with proton propagation 
limited by the GZK effect, but with neutrinos 
able to propagate over cosmological distances. 
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